arXiv:hep-lat/0112014vl 11 Dec 2001 


The extraction of qa from finite volume systems: The long and short of it 


Thomas D. Cohen 

Department of Physics, University of Maryland, College Park 20748 
UMPP#02-085 DOE/ER/40762-247 

Due to a pion-pole contribution, the nucleon axial current matrix elements can be visualized in 
position space as becoming delocalized as the chiral limit is approached, with one-third of the current 
at distance of order l/m,r- However, this delocalization effect will not cause calculations of qa in a 
finite system with standard boundary conditions (periodic for boson fields, antiperiodic for fermions) 
to have large finite volume effects. This is seen by calculating qa using a nonstandard current with 
axial quantum numbers. The matrix elements of this new current lacks a pion pole and is not 
delocalized; however, its zero-momentum matrix element is identical to that of the standard axial 
current both for infinite and finite volumes. 


The coupling of the nucleon to axial currents at zero 
momentum transfer is hxed by the coupling constant qa- 
Due to chiral symmetry and the key role played by pions 
in the low-lying dynamics of the nucleon, understanding 
how QA emerges from the underlying dynamics of QCD is 
essential if one wishes to understand how the low energy 
dynamics of the nucleon is connected to QCD. Attempts 
to calculate gA in lattice QCD have given values which 
are well below the experimental value Q . It has recently 
been argued that this failure may be due to very large 
finite volume effects [^. One may ask whether there is 
a basis for understanding the emergence of large finite 
volume effects from the structure of QCD. This issue was 
brought into sharp focus in a recent paper by Jaffe [^. 
Jaffe made an elegant argument on very general grounds 
that the spatial distribution of the axial current in the 
nucleon is necessarily spread over large distance as the 
chiral limit is approached, a phenomenon he denotes by 
“delocalization”. In particular he showed that as the 
quark mass goes to zero, one-third of the axial current is 
pushed off to infinite distance. For finite quark mass, this 
implies that of order one-third of the axial current can be 
found at a distance of order m~^. It was suggested in Ref. 
1^ that this effect could lead to large hnite volume effects. 
Intuitively it is not implausible that this delocalization of 
the axial current could lead to large hnite volume effects 
in calculations of gA'- if the pion is light enough so that 
large amounts of the axial current are pushed to regions 
which are beyond the lattice volume one might expect 
that the lattice calculation will miss these contributions 
leading to substantial errors. 

The purpose of this paper is to observe that although 
ref 1^ is correct that one-third of the axial current is 
delocalized in the chiral limit, the intuition that this 
implies large hnite volume effects is not (provided the 
standard boundary conditions of periodicity for bosons 
and anti-periodicity for fermions are imposed). To make 
this statement a bit more precise it is necessary to look 
a bit more carefully at the heuristic argument why one 
might expect large hnite volume effects. The basic issue 
is the play off between the inhnite volume and chiral lim¬ 
its. The most serious worry is that as the chiral limit 


is approached, one-third of the axial current will be lost 
outside the box for any hnite sized box. Thus to get 
convergence of gA one needs to choose a box size signif¬ 
icantly larger than one with a length scale set by llm-n- 
Mathematically this would be rehected by nonuniform 
limits of the calculated value of gA as one approaches the 
inhnite spatial volume and chiral limits. Dehne to 
be the matrix element of the axial current between nu¬ 
cleon states as calculated with hxed quark mass (which 
hxes m,r) in a hnite spatial volume. The principal con¬ 
cern is that calculated in the limit R ^ 00 , mT^ 0 
will only be two thirds of calculated in the limit 
m-n ^ 0, R —> 00 . However as will be shown below, if 
the standard boundary conditions are imposed the limit 
should be uniform: 
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Thus, the major intuitive concern about large hnite vol¬ 
ume effects can be put to rest. There is a secondary 
worry, however, that if much of the strength of the un¬ 
derlying integrals come from regions comparable to the 
box size, the dynamical effects induced by the bound¬ 
ary could alter this part of the distribution even if the 
strength is not actually lost. This could lead to hnite vol¬ 
ume effects which are larger than those for typical baryon 
matrix elements. This secondary concern has not been 
stated in precise form and, hence, it is hard to precisely 
contradict. However, it will also be shown below there 
is no reason based on the delocalization of current dis¬ 
cussed in Ref. js) to expect that hnite volume effects to 
be larger for gA than those for typical matrix elements. 

Let us begin with a brief discussion of the demonstra¬ 
tion that the axial current is delocalized. In a certain 
sense, the argument for the delocalization is quite old; 
the result is well known in the context of mean held chi¬ 
ral soliton and bag models of the nucleon [Q. In these 
models the nucleon is treated as a localized source and 
the pion held is treated classically. The localization of 
the source and the mean held treatments are typically 
justihed from large Nc QCD arguments. The axial cou¬ 
pling constant is calculated by integrating the axial cur¬ 
rent over the nucleon and then matching to the general 
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Lorentz covariant expression for axial matrix elements 
assuming that the nucleon source is localized. The result 
is that as TO,r 0 the integrals which give gA receives 
one-third of its contribution from distances of order tUtt, 
and as one approaches the chiral limit one-third of the in¬ 
tegrated strength is pushed off to infinity. In a sense, the 
crux of the recent argument of Ref. is that these well- 
known results in the context of mean held chiral models 
can be translated into a general model-independent re¬ 
sult with essentially no change in structure. 

To begin, consider the general structure of axial matrix 
elements of the nucleon in terms of form factors: 

(fV(p^s',J')K(0)|IV(p,s,J3))= (2) 

+ hA{q^)q^lb) u(p,s,h) , 

where u(p, 3 , 73 ) is a free isospin one-half Dirac spinor 
with a third component of spin, s, and a third compo¬ 
nent of isospin, I 3 |^. The form factors, gA{q^) and 
hA(q^), are generally referred to as axial and induced 
pseudoscalar form factors, respectively. The axial cou¬ 
pling constant, gA, is dehned as the axial form factor at 
zero momentum transfer. The key point is that the axial 
current acting on the vacuum can create a pion which 
implies that the induced pseudoscalar form factor has a 
pion pole contribution: 

hA{q^) = h^:^^{q^) + hf°'^\q^) 

= -2f^g^NN , (3) 

where is the pion decay constant (with the convention 
that /tt « 93MeV), g.„NN is the pion-nucleon coupling 
constant, and the non-pion-pole part of the in¬ 

duced pseudoscalar form factor. This non-pion-pole con¬ 
tribution to Ha is denoted as “short” to indicate the scale 
over which it varies substantially is the hadronic scale, 
A ~ IGeV, which is much larger than TOtt and hence 
corresponds to a spatial extent which is much smaller 
than The pion-pole contribution to the induced 

pseudoscalar coupling, however, varies substantially if 
changes on a scale comparable to which implies con¬ 
tributions on a spatial extent of order This rapid 

variation of the pion pole term is the origin of the effec¬ 
tive delocalization of the axial current. 

Rather than recapitulate the demonstration of Ref. ||] 
in detail, an equivalent and somewhat simpler derivation 
of the same result will be presented here. Consider a spa¬ 
tial component of the axial current between spinors. (A 
spatial component is used since the temporal component 
vanishes at zero momentum transfer). At hnite momen¬ 
tum transfer the results are frame dependent. It is natu¬ 
ral to work in the Breit frame of {piniuai = -pfinal)- The 
Breit frame has the virtue of having momentum transfer 
but no energy transfer. Thus the Fourier transform of 
the matrix element in the Breit frame may be directly 
interpreted in terms of position space; the relative time 


variable plays no role. Introducing the symbol G{q) as 
the matrix element of the third spin and isospin compo¬ 
nents of the axial current between spin-up protons in the 
Breit frame and a three-momentum transfer of Q, and 
using Eq. (^) and standard properties of Dirac spinors 
one has: 

G(Q2) = {n{Q, 1/2, l/2)|Ai(0)|iV(-Q, 1/2,1/2)) = 

[gAi-\Q?) (l + ^6^)+/^A(-|QnS) 

VQ2-H4M2+2M 
^ 4M ’ 

where M is the nucleon mass. The choice of the third spin 
and isospin components and the proton spin up states 
was arbitrary and done for the purposes of specificity; 
other choices of components and spin and isospin projec¬ 
tions of the state lead to identical results. 

From Eq. (^) it is immediately apparent that 
= 9 A, providing that there is no singu¬ 
larity in the induced pseudoscalar form factor at = 0. 
For any finite TOtt this is true and the limit both exists 
and is uniform in the sense that one can take the Q —> 0 
limit from any direction and will get the same answer. 
Because the limit Q —> 0 limit is uniform, if we define an 
angular average of G{Q), 

G(Q2) = ±J dflG(Q) , (5) 

then it is immediately clear that for any nonzero m.^-, 

= gA ■ ( 6 ) 

We wish to study G{Q^) for both g and m.,^ well below 
the typical hadronic scale of A ~ I GeV. Assuming 
gA{q^) and hA{q^) both vary over scales of order A, one 
sees from Eqs. (2), d), d) and (I) that 

G{Q')=9a - +0«3VA“,mJ/A“) . 

AM G -I- TO* 

(7) 

Imposing conservation of axial charge in the chiral 
limit on Eq. (^) yields the Goldberger-Treiman relation 
9 aM = fTr9-!rNN + Oiml/A^) ® so that 

G{q^) = gA fl- —2-) +C)(gVA^TO^/A2) . (8) 

V 3q^ + mi J 

Note that the preceding equation is the leading order 
term in chiral expansion; while there are correction terms 
of order Q^/A^ and toJ/A^, mi has been included to 
all orders. 

The delocalization of axial charge discussed in Ref. 
is clear from Eq. (||). The term proportional to G^/(Q^-I- 
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TO^) clearly contributes 1/3 of the (5 = 0 result. It 
is also clear that this term varies rapidly over momen¬ 
tum scales corresponding to spatial variations over long 
distance scales ~ l/m^. As ^ 0 the spatial dis¬ 
tance over which this term contributes becomes arbitrar¬ 
ily large, i.e. the current becomes delocalized. To under¬ 
stand this as a spatial delocalization of the current in the 
chiral limit, one needs a method to visualize the current 
in position space. The natural way is to define a spatial 
axial current density pa{t) as the Fourier transform of 
G{q^). If the current is delocalized in the chiral limit, 
then all of the positive moments of pa{t) will diverge. 
These moments are given by 

{r'‘")A = ^jd‘rr^yAr) 

= (s) 

Using Eq. one has 


whose eigenstates have a well-defined energy. These en¬ 
ergy eigenstates are most easily treated as being box- 
normalized. Because these states are associated with 
fields which have well-defined properties under trans¬ 
lation in the x, y, or z direction by L, the energy 
states are all eigenstates of the momentum operator with 
p = 2TrL~^{nxX + riyp + n^z) (integer n’s) for bosons 
and with p = 'KL~^{{2nx + I)x -I- {2ny + l)y + {2nz + 
l)z (integer n’s) for fermions. The box breaks rota¬ 
tional symmetry and hence angular momentum is no 
longer a good quantum number. However, the spin de¬ 
gree of freedom is still well defined for fermions in the 
sense that for every momentum there are two degener¬ 
ate eigenstates. Thus the nucleon states can be labeled 
as Ip,/ 3 ,s) with the normalization (p',/'a, s'|p,/s, s) = 
Sl>3j3Se>,sSn\,n^Sn’y,nySn'^,n^- One can define the finite 
volume version of pA as the diagonal matrix element of 
the third space and isospace components of axial current 
between spin up proton states on minimum momentum 
appropriately normalized to the volume of the box: 


{r^^A = 


5An!n,=i,„(2 + 4/) 






for n > I . 


( 10 ) 


pVoi = L^{pJ:i^s\Al{x,t)\pj3,s) 

= {p,h,s\f d^xAl{x,t)\p,l3 ,s), (II) 

♦/box 


all of which diverges ~ as the limit ^ 0 is 

approached. This indicates a delocalized current in the 
chiral limit. 

Thus, one-third of the axial current density can be vi¬ 
sualized as being delocalized as the chiral limit is ap¬ 
proached. While this represents an important insight into 
the physics, it is important to recall that the underlying 
dynamics depends on nothing beyond the fact that the 
axial current matrix element has a pion pole contribu¬ 
tion. In this paper, the central issue is the relationship 
between the delocalization of the axial current and ef¬ 
fects due to finite spatial volumes. Of course, real lattice 
calculations have a number of sources of errors including 
statistical and discretization errors as well as errors due 
to finite spatial and temporal extent. Here, the focus is 
on finite spatial volume effects. To isolate these effects 
it is useful to consider the behavior of continuum QCD 
on a finite three-dimensional volume with infinite tem¬ 
poral extent. From the prospect of lattice QCD this can 
be thought of as working with a sufficiently fine mesh, a 
sufficiently large number of lattice points in the time di¬ 
rection and a sufficiently large number of configurations 
so that discretization errors—errors associated with finite 
temporal extent—and statistical errors are all negligible. 

The calculation of pA in finite spatial volume QCD 
is straightforward in principle. We consider the three- 
dimensional space to be a cube of length L, and im¬ 
pose the boundary conditions that fermion fields are 
anti-periodic and boson fields are periodic. Since the 
time direction is infinite, one can define time evolution 
as being generated by a hermitian hamiltonian operator 


where the second form of Eq. ( 0 ) stems from transla¬ 
tional invariance of the matrix element with these bound¬ 
ary conditions. 

Before proceeding it is useful to distinguish between 
the method of calculating pA for the finite spatial sys¬ 
tem discussed above as an integrated quantity in a state 
of good energy from the method actually used in lattice 
calculations. In lattice calculations one never calculates 
exact states of good energy. Nucleon matrix elements 
such as PA are calculated as the ratio of a three-point cor¬ 
relation function to a two-point function with the times 
of the three currents well seperated in time. As the seper- 
ation in time between these currents becomes large, the 
system becomes dominated exponentially by the lowest 
energy state with quantum numbers which overlap those 
of the initial and final currents. As the time seperation 
between all three currents goes to infinity the ratio of the 
three-point function to the two-point function asymptoti¬ 
cally approaches the matrix element in the nucleon state 
considered above. Thus the matrix derived using the 
states of good energy are equivalent to the ratio of the 
three-point functions to the two-point functions used in 
lattice calculations proven that the times are arbitrarily 
long. If the times are not arbitraily long there are finite 
corrections. 

The issue of concern is whether the delocalization of 
the axial current for the infinite volume theory leads to 
large finite volume errors. It is easy to show that this 
is not the case by considering the matrix element of a 
new axial current operator. Suppose that we can find a 
current operator which in the infinite volume limit does 
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not suffer from the delocalization effect seen in the con¬ 
ventional axial current. If it can be shown that the in¬ 
tegrated matrix element of this new current in the finite 
volume theory precisely equals g\°^, then one can replace 
the original finite volume calculation of qa with one based 
on the new current. Doing this has the virtue that de- 
localization does not occur in the new calculation and 
hence cannot be responsible for large finite volume ef¬ 
fects. Since equals the finite volume QCD matrix 
element of the new current integrated over space, it fol¬ 
lows that also cannot have a large finite volume effect 
due to delocalization. It is straightforward to find a cur¬ 
rent that has these properties. Since the delocalization 
is traceable to the pion pole, we need to find a current 
whose nucleon matrix elements lack such a pole. 

Consider the current A'^ defined by 

+ , 12 , 

where is the D’Alembertian operator. Although 

this current has no well-defined chiral limit, it is well 
defined for any finite value of the quark mass and this is 
sufficient for the present purpose. The matrix elements 
of this current can easily be expressed in terms of the old 
one: 

IV(p',s',/')|A'^(0)|iV(ps,/3)) = 

u{p, s', l3)^{Ta {g' + h' A[q^)q'"lb) u{P: S, I 3 ) 

5 A( 9 ^) = (1 - 9Aiq^) 

h'A{q^) = (1 - hA{q^) ■ (13) 

It is clear from this form that as g —> 0, g 'a gA and 
h'A —^ 0. For the infinite volume case, one can define the 
function G"(Q^) in direct analogy to G{Q'^), 

= ^ J dDG'(Q) 

G'iQ) = {N{Q, 1/2, l/2)|A'^(0)|iV(-Q, 1/2,1/2)) . (14) 

Inserting the matrix elements for A' from Eq. di) 
into the definition for G'{0‘^) gives G'{q^) = (1-1- 
/m‘^)G{Q^). Since, Eq. (^ gives G(g^) to lowest chi¬ 
ral order {i.e. up to corrections of order rn^/A'^ or j Is? 
but including Q'^Ini^ to all orders) we immediately have 
the lowest order expression for G'{Q'^)-. 

^(Q2) = sa (1 + + G(QVA^ ml/A^) . (15) 

^From the form of Eq. a it is immediately clear that 
for any finite quark mass, G'(0) = G(0) = gA- Thus, in 
principle one can use either A'“ or A“ to extract gA- 
The form of Eq. ( IT^ shows that the current associ¬ 
ated with A'/ is not delocalized in the sense that a finite 


fraction of the integrated strength is at distance scales 
of order l/m,r- Were it so delocalized, all positive mo¬ 
ments of the distribution would diverge as ^ 0 with 
~ 1/mJ" . These moments are given by 

== n (2 +• <>«) 

which from Eq. ( p^ ) implies that 
(?’^)a' = 2 

< G(m-2") for n > 1 . (17) 

Thus none of the higher moments diverge with in a 
fashion consistent with a delocalized distribution. 

One might ask how the moment could diverge with 
1 /toJ if the distribution is localized. The answer is sim¬ 
ple: the spatial distribution pA'{r) could be of the fol¬ 
lowing form: 

A® 

PA'{r) = gAA^fi{Ar) H-(18) 

ml 

with J dr 47 rr^/i(r) = 1 

where /2 (x) goes to zero as x ^ 00 faster than any power 
law. The first term in Eq. (^^ integrates to the full value 
of gA while the second term integrates to zero due to the 
derivative. Moreover, the second term, although of large 
overall scale, is of very short range. Thus it leads to 
the moments of Eq. ( 0 )- It is clear that as the chiral 
limit is approached the second term remains sizable only 
over very short distances and cannot be associated with 
delocalization. 

On the other hand, it is straightforward to see that if 
one uses A'“ to calculate gA on a finite volume one gets 
exactly the same result as if one uses A“. To see this 
consider 

9a°' = {Pl ^' 3 , s'l /|^^^d3xA'3(x, t)\p, h, s) 

= (p'./'a.s'l/i^^^d^x + Al{x,t))\p,h,s) . (19) 

The second equality shows that the matrix element of 
A'\ has two parts—the matrix element of Ag and a total 
derivative term. This second term can be written as 

(p',/' 3 ,s'| / d^x ^^—5 -A|(x,t = 0 ))|p,/ 3 ,s) (20) 

^boX ^TT 

= —{p'J'z,s'\ / A^xn-\IAl{x,t = Q)\p,h,s), 

^TT X surface 

where n is an outward oriented normal vector; the time 
derivative vanishes since our matrix element is evaluated 
for an energy eigenstate and the divergence theorem is 
used to replace the volume integral by a surface integral. 
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On the other hand, the standard boundary conditions 
implies that A\{x,t) = A^{x -\- Lx A) = A\{x + Ly) = 
A\{x + Ly) which in turn means that VA\{xA = 0) on 
one wall of the box is identical to VA\{xA = 0) on the 
opposite wall. Since the outward normal is opposite on 
opposite walls the integral vanishes. Only the first term 
of Eq. (|^) contributes and thus one has 

This completes the demonstration. The current 
unlike remains localized in the chiral limit and hence 

H- 

cannot suffer from large finite volume effects due to de- 
localization. However, it gives the same result for gA on 
a finite system as the original current provided the stan¬ 
dard boundary conditions are imposed. Thus gA calcu¬ 
lated on a finite volume from also cannot have large 
finite volume effects due to delocalization. The proceed¬ 
ing demonstration shows that the intuition that the delo¬ 
calization of the axial current implies large finite volume 
effects is not correct. 

It is useful to understand in a more heuristic manner 
why this intuition fails. Let us recall that the underlying 
reason for concluding the current was delocalized in the 
chiral limit was the rapid variation of G'(Q^) for small 
pion masses assoicated with the pion pole. To see what 
is happening, it is sensible to restate the intuition that 
large finite volume effects might be expected in terms of 
momentum space. If one visualizes restricting the nu¬ 
cleon to a box, one might think that the states will not 
have well-defined momenta and the characteristic mo¬ 
menta in the state will be of order 1/L. In such a case 
the matrix element would have components of momen¬ 
tum with Q ^ 1/L and the rapid variation of G{Q^) 
might be expected to yield large finite volume effects un¬ 
less m.n.L >> 1. However, this intuition is faulty. In par¬ 
ticular, the perception that the momentum is not well- 
defined, while valid for boundary conditions in which the 
fields vanish on the walls, is not valid with the standard 
boundary conditions (periodic for bosons; antiperiodic 
for fermions). As discussed above, however, the stan¬ 
dard boundary conditions are compatible with states of 
good momentum. The effect of imposing the standard 
boundary conditions does not mix momenta. Instead 
these boundary conditions allow momentum eigenstates 
but restrict the allowable momenta. Thus, with these 
boundary conditions one can have Q = 0 exactly, and 
the fact that matrix elements vary rapidly with Q near 
(5 = 0 becomes irrelevant. The heuristic argument why 
delocalization in the chiral limit does not lead to large 
finite volume effects depends critically on the use of the 
standard boundary conditions. This is consistent with 
the more formal demonstration above in which standard 
boundary conditions were also essential in showing that 
the surface integral vanishes. 

It is also useful to understand the failure of the in¬ 
tuition that for finite volumes and a delocalized source 
one expects to lose some of the strength out of the box 
and hence to have large finite volume effects. Again the 


boundary conditions provide the explanation. The key 
point is that with these boundary conditions, a nucleon 
in the box is equivalent to an infinite system with a cu¬ 
bic array (with a lattice spacing of L ) of boxes each 
containing a single nucleon. Thus, while strength can be 
lost from a box into its neighbors due to delocalization, 
a compensating strength may be gained from strength 
escaping into the box from its neighbors. 

In summary, it has been shown that the delocalization 
of the axial current discussed in Ref. does not lead to 
large finite volume effects. This was shown explicitly by 
constructing a current which is manifestly not delocalized 
and which leads to the same finite volume value for gA 
as the standard current. The demonstration depended 
on the imposition of the standard boundary conditions. 
It was argued at a heuristic level why it is not surpris¬ 
ing that delocalization does not yield large finite volume 
effects provided these boundary conditions are imposed. 

This result leaves open the question of why present 
lattice calculations substantially underpredict gA- Re¬ 
cent lattice results do suggest large finite volume effects 
[^. This does not contradict the present result. Here 
it is only shown that the delocalization effect will not 
lead to large finite volume corrections. Of course, other 
effects can cause big finite volume effects. Given the rel¬ 
atively small volumes used for the present generation of 
lattice calculations it is not too surprising that large fi¬ 
nite volume effects are present. However, there is also no 
reason to connect them with the delocalization effect of 
Ref. [|j. Ultimately, reliable lattice calculations should 
become available on larger lattices and the playoff of the 
chiral and infinite volume limits can be explored in the 
calculation of gA- iFrom the arguments presented here, 
one expects such calculations to indicate that the com¬ 
bined limit is uniform. 
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